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摘 要 考虑 营 发 冷却 效应 、 气 相 物 质 扩散 、 气 液 两 相传 热 及 界面 变形 , 采用 任意 拉 格 朗 日 - 欧 拉 法 数值 模拟 研究 了 接触 
fi 0. RRR Ec. Ma 数 对 0.65cSt 硅油 液 滴 固定 接触 角 模 型 蒸发 及 Marangoni 对 流 的 影响 . 发 现 : 液 滴 无 量 纲 体 积 V* 
变化 与 接触 角 09、Ma 数 有 关 , 而 与 Bc 数 无 关 ; H V*2/3 随时 间 线 性 减少 . 当 9 < 0.175 rad 时 , 随 着 Ma 数 的 升 高 , 液 
滴 内 部 依次 出 现 热 毛 细 涡 流 结构 、 单 个 BM 涡流 结构 及 多 个 BM 涡流 结构 . 当 0.175 rad< 0 < 0.479 rad 时 , 存在 热 毛细 
涡流 结构 或 多 个 BM 涡流 结构 或 两 种 同时 存在 ,而 无 单个 BM 涡流 结构 . 当 0 > 0.479 rad 时 , 仅 存 在 热 毛细 涡流 结构 . 
当 0.105 rad< 0 < 0.401 rad AY, 每 两 种 涡流 结构 的 分 界 Ma 数 均 随 9 增 大 而 减少 . 当 0 > 0.401 rad 时 ， 热 毛细 涡流 结 
构 与 多 个 BM 涡流 结构 的 分 界 Ma 数 随 9 增 大 而 增 大 . 
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= Abstract Taking into ace bf evaporation cooling, mass diffusion, two-phase heat transfer and 
© gas-liquid interface aie 


contact angle on a heat 


ion, the evaporation of a droplet of 0.65cSt silicone oil with constant 
substrate is numerically investigated using AWbitrary Lagrangian-Eulerian 
method. The effects of contact angle 0, evaporation number Ec ad iarann number Ma on the 
evaporation and the Marangoni convection are analyzed. N (ital results find that the variation 
of volume V* is affected by 6 and Ma but not by Ec, a e V*?/3 varies with t* linearly. With 
SA Bénard-Marangoni convection (BM 


increasing Ma, the thermocapillary convection, singl 

a sequence for 0 < 0.175 rad. Thermocapillary 
th of them occur for 0.175 rad< 0 < 0.479 rad 
but without single-vortex ee onvection. Only thermocapillary convection appears 


convection) and multiple-vortex BM convection 
b 


convection or multiple-vortex BM convection 


for 0 > 0.479 rad. The neutral Ma between two kinds of vortex convection decreases with increasing 
contact angle for 0.105 rad< 0 < 0.401 rad, while it increases for 6 > 0.401 rad. 
Key words evaporating droplet; two-phase heat transfer; interface deformation; Marangoni convec- 


tion 


0 引 F 


液 滴 落 发 在 工程 实际 中 的 运用 十 分 普遍 , 例如 : 
相 变 冷却 技术 出 、 喷 雾 薄膜 涂 层 工艺 四、 打印 着 
色 工 艺 回 、DNA 芯片 制造 和 四、 DNA 分 子 调 整 及 拉 
伸 本、 血液 沉积 疾病 诊断 技术 中 等 . 由 于 液 滴 为 弧 
形 表面 , 导致 液 滴 表面 总 存在 切 向 和 法 向 温度 梯度 . 
切 向 温度 梯度 引起 热 毛细 对 流 ， 足 够 大 的 法 向 温度 
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梯度 引起 Bénard-Marangoni 对 流 ( 即 BM 对 流 ), 文 
献 中 不 加 区 分 时 两 者 均 称 为 Marangoni 对 流 。 

液 滴 蒸发 及 其 内 部 Marangoni 对 流 的 研究 大 多 
集中 于 固定 接触 线 模型 (CR 模型 ) !'- 引 , 包括 液 滴 蒸 
发 过 程 的 质量 变化 (蒸发 率 )、 几 何 形变 (接触 角 、 液 
滴 高 度 、 体积 )、 沉积 现象 、 表 面 温度 以 及 Marangoni 
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对 流 及 其 失 稳 等 . 其 中 文献 [9] 采用 归 一 化 法 研究 了 
几何 形变 、 液 滴 寿 命 与 润 湿 半 径 、 基 底 温 度 、 环 境 
湿度 的 关系 。 实际 上 ， 由 于 基底 表面 特性 不 同 ， 液 
WAAR ST CR 模型 外 ， 还 有 接触 角 保 持 不 变 的 情 
况 OO), 即 固定 接触 角 模型 (CA 模型 ), 而 目前 对 CA 
模型 下 液 滴 蒸 发 及 其 内 部 Marangoni 对 流 的 研究 相 
对 较 少 。 已 有 研究 主要 采用 理论 和 实验 的 方法 ， 研 
究 了 液 滴 质量 Lt1-1 引 、 液 滴 体 积 OS) 变化 规律 及 热流 
体 波 ( 失 稳 的 热 毛细 对 流 ) 的 波 数 与 初始 液 滴 几 何 尺 
寸 、 时 间 的 关系 M, 在 “one-side” 模 型 下 , 文献 [15] 
给 出 了 Bénard-Marangoni 对 流 模拟 结果 , 但 并 没有 
深入 研究 ， 

为 些 ， 本 文采 用 任意 拉 格 朗 日 - 欧 拉 法 ,数值 
模拟 固定 接触 角 模 型 情况 下 加 热 平 面 上 营 发 液 滴 内 
Marangoni 对 流 ， 主 要 分 析 接 触角 0. Ec HL. Ma 数 
对 它 的 影响 。 研究 对 理解 固定 接触 角 营 发 液 滴 内 部 
流动 和 传 热 特征 具有 指导 意义 。 


1 数学 物理 模型 


1.1 物理 模型 He 
物理 模型 如 图 1 所 示 ， ne eGR, 
初始 润 湿 半 径 Re = 2.5 mm, 0.105 ~ 


接 角 

0.524 rad; 基底 长 度 为 1.5R。， pen 65Re. 

MWS BRR ABM, Re ae 散 、 气 液 两 相 流 

ae he 两 相 工 质 分 别 为 空 
气 、0.65cSt 硅油 ,其 物性 参数 如 表 1 Arm. 


Axis 


0.65 cSt silicone oil 


图 1 基底 上 液 滴 蒸发 物理 模型 


Fig. 1 Schematic model of evaporating liquid droplet on a flat 


substrate 
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表 1 25 °C 时 0.65cSt 硅油 、 空气 


的 物性 参数 【9 


Table 1 Physical properties of 0.65cSt silicone oil 
and air at 25 °C 06] 


2 


参数 0.65cSt 硅油 空气 

密度 p/(kg/m?) 760 1.18 
BARRE u/((N-s)/m?) 4.94x10-4 1.83x10-5 
导热 系数 和 /(W/(m:K)) 0.1 0.01 
比热容 cp/(J/(kg:K)) 2000 1004 
表面 张力 系数 o/(N/m) 1.52x10-2 一 
表面 张力 温度 系数 oT/(N/(m:K)) 9x1075 = 
RRM Ht/(J/kg) 2.23x105 

浓度 扩散 系数 D/(m?.s) 5.60x10-6 -— 


1.2 数学 模型 及 计算 方法 
模型 采用 轴 对 称 模型 , 基本 假设 包括 : 液 滴 初 始 

形状 为 球 冠 形 ; 除 表 面 张 力 随 温 度 变化 以 外 , 其 余 物 

性 参数 均 为 常数 ; 基底 表面 为 等 温 边 界 条 件 ; 表面 张 

力 在 流 中 占 主导 作用 (Bo ~ 0.01)， 忽 略 重 力 

HIERS 钨 液 界 面 蒸 汽 浓 度 为 饱和 浓度 ; 气 、 液 两 相 
se 可 压缩 流体 . 考虑 到 VOF 法 难以 计算 变 

乡 较 大 的 问题 、 而 Level Set 法 不 满足 质量 守恒 ， 

文 借鉴 文献 [17] 的 计算 方法 , 采用 任意 拉 格 朗 日 - 欧 

拉 方 法 (ALE) 追踪 气 液 界面 变形 . 

1.2.1 无 量 纲 控 制 方程 


质量 守恒 方程 : 
V-u*=0 (1) 
动量 守 a 
a A 
Set (u*-V)u*=- Vp I + —V2u* (21) 
\\\ 
气相 
oe (uz - V) u*= — at +—V?u* (2-2) 
能 量 守恒 方程 
Bt (ux. V)T Po wh (3-1) 
气相 
OT* s 1 * 
Bt (uz . V) Per, VT (3-2) 
气相 扩散 方程 
Oc* ‘ el ‘ 
ae + (te Vic =p Ve (4) 
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其 中 ww 为 流体 流动 速度 与 网 格 移动 速度 之 差 : w: = 
u* — Umesh 
1.2.2 边界 条 件 
1) 流动 边界 条 件 
a) 对 称 轴 采 用 对 称 边界 条 件 
液 相 : 


u*-n=0, -7 十 five -n=0 (5-2) 
b) 气 液 界面 力 平衡 条 件 
法 向 方向 : 


1 
-ve :nn 一 Ler 十 mre : nn 二 
el 1 
g 


mir oe” 
-» | [Ke] m 
(t) vi (6-2) 


气 液 界面 质量 流量 与 气 液 相 速度 关系 I, 


(uig gs — Vlg_1) = (+ = 1) (n- Vc") 
Ps 
其 中 ,“: ”为 张 量 双 点 积 运算 符 ,“m” “g” 
面 外 法 向 单位 向 量 和 切 向 单位 向 量 , “nn”, “nt” 
它们 的 并 矢 。 曲率 k* = Re/re, re HERS. 
c) 液 固 、 气 固 界面 采用 Naver-slip 滑 移 边界 条 


件 
液 固 界面 : 
[pr + aye] : nt = TH (8-1) 
气 固 界面 : 
Ea + aye] : nt = Tee -t (8-2) 


其 中 , 6 一 DB/ Re 为 无 量 纲 滑 移 长 度 , 6 ARBRE. 


d) 气相 外 边界 采用 定 压 边 界 条 件 
Po =0 (9) 


2) 传 热 边 界 条 件 
a) 对 称 轴 采 用 对 称 边界 条 件 


—n-(VT"*) =0 (10) 
b) 气 液 界面 能 量 平衡 边界 条 件 


n. [OvT (V7) =Ec(—n-(—Vet)) (11) 


c) 气相 外 边界 采用 恒温 边界 条 件 
T* =0 (12) 
KE 温 边 界 条 件 
Te 一 1 (13) 


A 
a) 对 称 轴 采用 对 称 边界 条 件 


-n (Vč) =0 (14) 
b) AWF ERA WMR BE FR RF 
Chg = (15) 
c) 气相 外 边 缀 膝 用 恒定 浓度 边界 条 件 
. O c =0 (16) 
ONREG 
i -n (Vč) =0 (17) 


RG 1.2.3 初始 条 件 


Ut=0 = 0, Dt-0 = 0, Cho =0 (18) 


方程 中 几何 尺寸 、 速 度 、 时 间 、 压 力 、 密 度 、 
黏度 、 导 热 系 数 的 无 量 纲 尺 度 分 别 为 Re, we、 Re/ Ue, 
puž, pir m, 入， 其 中 we = 2rMD(csat — Coo) /piRe. 
无 量 纲 温 度 定义 为 7* = (T — Too) /(Tw — Too), 无 
量 纲 浓度 定义 为 * = (c 一 cwo)/(csat 一 cw); FER 
s、lg-l、lg-g、1、g、w、o0 分 别 表示 气 液 界面 、 界 面 
液 相 侧 、 界 面 气 相 侧 、 液 相 、 气 相 、 基 底 、 气 相 外 
边界 ; 浓度 Coat. Coo 满足 理想 气体 状态 方程 Csat = 
psat/(RT)、coo = Psat RH/(RT); 压力 Psat (单位 kPa) 
是 通过 安 托尼 方程 ng lg(psat) = A- B/(T + C) 得 
到 ; 其 中 A = 7.07406, B = 1657.46, C = —46.13, 


To = 0, 


ChinaXiv 合 作 期 刊 
38 卷 


4 期 


无 量 纲 特 征 参 数 定义 为 : 雷诺 数 Re = uR-/ 
(u/p), 普 朗 特 数 Pr = Hcn/ 和 ,， 传 热 佩 克 莱 数 Per = 
RePr = ueRe/(Mpoco)， 表 示 对 流传 热 与 扩散 传 热 
ZE, PAMEK Pep = weRe/D， 表 示 对 流传 
质 与 扩散 传 质 之 比 ， 毛 细 数 Ca = Hue/a， 马 兰 戈 
JER Ma = or(Tw 一 Too) Re/(mAi/picpa), RRR 
Ec = HM D(csat 一 Coo)/(Lw 一 Too), BARBRA 
热 与 扩散 传 热 之 比 。M 为 摩尔 质量 ，D 为 扩散 系 
BX, Ru = coo/csat 为 气相 外 边界 相对 浓度 , RAG 
体 常数 ，H' 为 气 化 潜 热 ,o 为 表面 张力 系数 ,oT 为 
表面 张力 温度 系数 . 


2 网 格 无 关 性 检验 


为 了 选择 合适 的 计算 网 格 , 选用 网 格 太 二 为 7.5 
um, 9 um, 10 pm, 12.5 um 四 种 网 格 进行 计算 . 其 
中 液 相 网 格 采用 上 述 四 种 网 格 尺寸 构成 的 均匀 网 格 ， 
气相 网 格 尺 寸 由 气 液 界面 处 网 格 尺 寸 向 气相 外 边界 


以 1.15 倍 的 增长 率 逐 渐 增 长 , 最 大 网 格 不 超过 LON 


hm。 计 算 180 s 后 液 滴 润 湿 半 径 及 表面 平均 温 
果 如 表 2 所 示 ， ts i sh 
虑 计算 时 间 和 精度 影响 , 后 续 计 算 选 取 尺 & m 的 
网 格 。 


表 2 网 格 无 关 性 验证 ， 润 ; 
Table 2 Mesh dependenty | r the droplet with 


oe 


c = 2.5 mm 


R. = 2.5 mm 
液 相 网 润 湿 半 7 相对 ”表面 平均 ”平均 温度 相 
R/um 径 7/nm RŽ% 温度 /*C 对 误差 % 
7.5 1858.60 29.44 一 
9 1864.66 0.33 29.44 0.0034 
10 1865.05 0.02 29.44 0.0034 
12.5 1869.18 0.22 29.44 0.014 


人 tr Ra om ES 
bee 


拟 文献 [18] 中 悬垂 半球 形 0.65cSt BERRIN 
度 变化 . 图 2 为 本 文 模拟 结果 与 文献 结果 比较 ， 最 
大 相对 误差 2.4% 出 现在 初始 阶段 误差 较 小 ,说明 
目前 采用 的 计算 是 正确 的 ， 
3 结果 与 分 析 
3.1 液 滴 节 发 过 程 及 其 影响 因素 

如 图 3 所 示 ， 实 线 为 模拟 的 液 滴 体积 变化 ， 虚 
线 为 拟 合 曲线 . WRR V 的 2/3 次 方 随时 间 线 
性 减 小 ,这 与 文献 [LL] ~ [13] 的 理论 与 实验 研究 结果 
趋势 一 致 . 且 接 触角 9 越 小 ， 其 无 量 纲 蒸发 春 命 越 


短 (如 图 , 当 V*3 = 0 时, 9 从 小 到 大 对 应 的 女 分 
别 为 0.50、0.66、0.82、0.97). 
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0 4 8 12 16 20 24 
t/s 


图 2 模拟 的 液 滴 顶 点 温度 与 文献 [18] 对 比 结果 
Fig. 2 Comparison of apex temperatures of droplet between 


1D resultl!8] and present simulation result 
aA, 


0=0.26 rad 


< 
eh 


2 
5 
N 
a 


0.0 A A N S 
0.0 G02 0.4 0.6 0.8 1.0 
N r 
提问 接 创 角 下 流水 发 过程 V*2/3 随时 间 t* 的 变化 
\¥ 3 Time evolution of droplet volume V*2/3 for different 
. angles (rad) 


通过 对 适用 于 CA, CR 模型 的 理论 公式 HH 
Me = 2xf(0)D(csat (Tig) 一 Coo)r 进行 分 析 ， 发 现 液 
滴 体 积 变化 受 三 方面 的 影响 : 1) 接触 角 o, feo [3 
BE 0 减 小 而 减 小 ; 2) 润 湿 半径 +,， 润 湿 半 径 7 BERR 
发 的 进行 逐渐 减 小 ;3) 浓度 coat (Tig) (Coo 不 变 时 )， 
安 托 尼 方 程 09 揭示 浓度 csat(Tis) 与 表面 平均 温度 
呈正 相关 关系 ， 而 表面 平均 温度 在 蒸发 过 程 中 变化 
相对 缓慢 .分析 这 些 因素 发 现 ， 相 同 蒸发 条 件 (Ec 
Bt. Ma 数 相 同 ) F, 液 滴 体 积 变化 主要 受 液 滴 几 何 
形变 (0, r) 的 影响 ， 

根据 1.2 节 中 无 量 纲 特征 数 的 定义 , Ec 数 随 Ry 
的 增 大 而 减 小 ( 即 : 外 边界 浓度 co 增加 ， 浓 度 差 
Csat (Tw) 一 Coo 减 小 ， 燕 发 冷却 效应 减弱 )，Meau 数 随 
Tw 增 大 而 增 大 . 模拟 中 通过 改变 Ry. Tw 来 研究 Ec 
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Rl. Ma 数 对 液 滴 无 量 纲 体积 变化 的 影响 , 发 现 两 者 
的 影响 (如 图 4) 与 固定 接触 线 模 型 有 一 定 的 差异 。 
在 CR 模型 F, Ec 数 、Ma 数 对 液 滴 无 量 纲 体积 
变化 无 影响 . 而 对 CA 模型 ， 虽 然 Bc 数 对 液 滴 无 量 
纲 体积 的 变化 无 明显 的 影响 , 但 Ma 越 大 ， 液 滴 无 
量 纲 体积 变化 越 慢 (无 量 纲 蒸发 寿命 越 长 ). 这 是 因 
HEE Ma BOK, WOM AHI TER ety) (图 
4 中 Ma 由 小 到 大 对 应 的 时 间 尺 度 分 别 为 : 597.04 
s, 422.09 s, 312.42 s, 238.68 s); 同时 液 滴 内 部 流动 、 
温度 分 布 也 随 之 发 生变 化 (如 图 6 所 示 ). 


1.0 上 
0.9 
0.8 
0.7 
~ 0.6 

Ma=34615.46 

0.5 È Ma=69230.77 

[ Ma=103846.15 

04 上 ”一 一 Mo=138461.54 
0.3 


0.0 0.2 0.3 x 
t SR 
图 4 9 = 0.524 rad 时 Ma XW V* 随时 间 t* 
see 
Fig. 4 The effects of Ma on the evolution of droplet volume 
V* with constant angle 0 = 0.524 rad 


3.2 液 滴 内 Marangoni 对 流 
如 图 5(a) Bras, 在 相同 的 Ma 数 下 ， 当 接触 角 


0 = 0.524 rad 时 , 液 滴 内 部 存在 单个 热 毛细 涡流 , K 
面 流 动 沿 三 相 接 触 线 向 顶点 流动 ， 表面 温度 随 流 动 
方向 逐渐 降低 。 当 接触 角 0 = 0.262 rad IY, 液 滴 内 
部 中 心 区 域 存在 两 个 BM 涡流 , 如 图 5(b) ras, 这 
两 个 BM 涡流 在 三 维 空间 上 组 成 一 个 BM Wi; 而 
表面 最 高 温度 一 直 在 三 相 接触 线 处 ， 因 此 在 液 滴 边 
缘 始 终 存在 由 三 相 接 触 线 向 液 滴 顶 点 流动 的 热 毛细 
涡流 。 热 毛细 涡流 与 BM 涡流 间 存 在 局 部 冷 点 ,而 
BM 涡流 之 间 存 在 局 部 热点 . 另外 , 在 Marangoni 对 
流 作 用 下 , 液 滴 内 部 等 温 线 发 生 了 弯曲 ; HEE BM 
涡流 结构 的 液 滴 其 等 温 线 弯曲 更 加 严重 , 说明 对 流 
更 强 。 

在 蒸发 前 期 ， 液 滴 内 部 流动 由 热 毛 细 对 流 逐 渐 
转变 为 BM 对 流 ,， 且 中 心 处 BM 涡流 强度 比 靠近 热 
毛细 对 流 M 涡流 更 强 、 涡 胞 尺寸 更 大 ， 如 图 6 
所 示 a = 34615.46， 直 到 t* = 1.67 x 107? 时 
} 部 热 毛细 涡流 与 BM 涡流 的 分 布 才 趋 于 相对 
Me, M Ma = 138461.54, t* = 3.35 x 1073 时 内 部 


DQ) 涡流 结构 已 相对 稳定 . 因此 Ma 数 越 大 ， 液 滴 内 热 


毛细 涡流 与 BM 涡流 达到 相对 稳定 的 时 间 越 短 ;， 且 
中 心 处 的 BM 对 流 占据 液 滴 内 部 的 空间 越 大 . 

除了 蒸发 过 程 的 影响 外 ，Ma 数 也 会 对 液 滴 内 
部 涡流 结构 产生 影响 , 如 图 7 所 示 . Ma = 5.54 x 10° 
时 , 液 滴 内 部 形成 单个 道 时 针 流 动 的 热 毛 细 涡 流 结 
构 ; Ma = 9.00 x 10° 时 液 滴 内 部 出 现 单个 BM 涡 
流 ， 其 流动 方向 3 浆 毛 细 涡 流 流动 方向 相反 ( 顺 时 
FLATT), FET PA AT BM 涡流 结构 形成 一 
个 BM 涡 胞 ，Ma = 1.73 x 104 时 , 液 滴 中 心 区 出 现 
PRBN 涡流 结构 ,两 者 形成 一 个 BM 涡 胞 结构 , E 


AY 


xv 


Thermocapillary 

0.2 convection (a) r *=0.118 1.00 
A 0.99 
6 . i 0.97 
i 1 1 | Er: I 1 f TONN ES EN EE Cer ey CS AE Us Ce Cn el We Cnet ee et Ce 0.96 
1.0 0.8 0.6 0.4 0.2 0.0 0.2 0.4 0.6 0.8 1.0 0.95 
0.2 Benard-Marangoni 0.94 
Thermocapillary convection (b) r *=0.071 0.93 
N 0.1 convection ano 
0.0 0.91 
0.90 

1.0 0.8 0.6 0.4 0.2 0.0 0.2 0.4 0.6 0.8 1.0 

re 


图 5 Ma = 69230.77 时 液 滴 内 部 等 温 线 及 内 部 流动 情况 (a) 9 = 0.524 rad, (b) 0 = 0.262 rad 
Fig. 5 The temperature field, streamlines and isotherm temperature lines for Ma =69230.77. (a) 0=0.524 rad, (b) 0=0.262 rad 
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Fig. 6 The evolution of temperature field and streamlines 
(b) t* = 3.35 x 1073, (c) t* 31.67 


N) 
ae (b) 
| Ee 


1.0 0.0 


0.1 


0.0 


0. 


201705.00924v1 


chinaXiv 


0 0.2 0.4 wT 0.8 
ZAS 
R (c 
0.0 0.2 0.4 0.6 0.8 


r* 


» i o 


0.2 0.4 0.6 
1.0 0.0 2 4 ~ 06 
Š: 
A 


1.00 
0.95 
0.90 
0.85 
0.80 
0.75 


图 7 接触 角 0 = 0.157 rad 时 液 注 内 部 温度 分 布 及 内 部 流动 情怀 Y Ma = 5.54 x 103, (b) Ma = 9.00 x 103, 


(c) Ma = 1.73 x 104, 
Fig. 7 The evolution of temperature field, streamlin 
(b) Ma = 9.00 x 10, ( Se 


在 液 滴 腰部 形成 新 的 小 涡流 。 Ma = 2.08 x > AY, 
液 滴 内 部 出 现 三 个 BM 涡流 ， 其 中 靠近 热 毛 细 对 流 
的 BM 涡流 的 流动 方向 为 顺 时 针 , 而 BM 涡流 之 间 
相互 耦合 流动 。 

由 于 液 滴 为 弯曲 界面 , 液 滴 内 必然 同时 存在 法 
向 温度 梯度 和 切 向 温度 梯度 。 只 要 存在 切 向 温度 梯 
E, 就 会 产生 热 毛细 对 流 。 所 以 ， 当 Ma 数 较 小 时 ， 
液 滴 内 以 热 毛细 对 流 为 主 . 当 Ma 数 足够 大 , 对 应 的 
法 向 温度 梯度 超过 BM 对 流失 稳 的 临界 条 件 时 , BM 
对 流 就 会 产生 . 相 比 热 毛细 对 流 , BM 对 流 走 过 的 路 
径 较 短 , 传 热 效率 高 , 因此 , 随 着 Ma 数 的 进一步 增 
K, BM 对 流 逐 渐 占据 主导 , 表现 为 BM 对 流 涡 胞 数 
逐渐 增多 , 而 热 毛细 对 流 所 占 的 区 域 相应 缩小 ， 


ww 


constant angel 0 = 0.157 rad. (a) Ma = 5.54 x 103, 
1.73 x 104, ( 


d) Ma = 2.08 x 104 


为 了 分 析 接 触角 对 液 滴 内 Marangoni 对 流 稳定 
性 的 影响 ， 本 文通 过 一 系列 的 计算 , 确定 了 不 同 接 
触角 下 Marangoni 对 流 结 构 之 间 的 分 界 Ma 数 , 如 
图 8 所 示 。 

接触 角 介 于 0.105~0.524 rad 之 间 , 存在 三 种 不 
同 的 涡流 结构 , 即 : 热 毛 细 涡 流 结构 、 单 个 BM 涡流 
结构 及 多 个 BM 涡流 结构 。 当 接触 角 为 0.105~0.175 
rad 时 ， 在 同一 接触 角 下 ， 随 着 Ma 数 的 升 高 ， 液 
滴 内 部 依次 出 现 热 毛细 涡流 结构 、 单 个 BM 涡流 结 
构 及 多 个 BM 涡流 结构 三 种 对 流 结 构 ， 当 接触 角 为 
0.175~0.279 rad 时 , 液 滴 内 部 仅 存在 多 个 BM 涡流 
结构 ; 当 接 触角 大 于 0.279~0.401 rad 时 , 液 滴 内 部 
依次 出 现 多 个 BM 涡流 、 热 毛细 涡流 及 多 个 BM W 
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流 两 种 对 流 结构 , 而 不 会 产生 单个 BM 涡流 结构 ; 当 
接触 角 为 0.401~ 0.479 rad 时 , 在 同一 接触 角 下 , 随 
着 Ma 数 的 升 高 , 液 滴 内 部 依次 出 现 热 毛细 涡流 、 多 
个 BM 涡流 两 种 对 流 结构 ; 当 接 触角 大 于 0.479 rad 
时 ， 液 滴 内 部 仅 存在 热 毛细 涡流 对 流 结构 ， 当 接触 
和 角 为 0.175~0.401 rad 时 , 即使 在 很 小 Ma 数 下 液 滴 
内 部 也 会 产生 多 个 BM 涡流 结构 ， 因 此 蒸发 冷却 效 
应 可 以 诱导 BM 涡流 结构 的 产生 ， 

图 8 中 黑色 线 为 单个 热 毛 细 涡 流 结 构 与 单个 
BM 涡流 结构 的 分 界线 ,红色 线 为 单个 BM 涡流 结 
构 与 多 个 BM 涡流 结构 的 分 界线 ， 蓝 色 线 与 绿色 线 
为 热 毛细 涡流 对 流 结构 与 多 个 BM 涡流 结构 的 分 界 
B, 黑色 线 与 红色 线 的 分 界 Ma 数 随 接触 角 增 大 逐 
渐 减 小 , 这 是 因为 随 着 接触 角 的 增 大 , 液 滴 的 厚度 增 
加 , 液 滴 内 流体 受到 基板 的 黏 滞 力 相对 减 小 , 同时 液 
滴 内 部 导热 减弱 而 液 滴 燕 发 冷却 效应 增强 ， 即 使 在 
相同 的 基板 温度 下 , 液 滴 内 Marangoni 对 流 也 更 强 . 
因此 ， 随 着 接触 角 的 增 大 ， 分界 Ma 数 逐 渐 减 小 . 
类 似 地 ， 蓝 色 线 与 0.279~0.401 rad 之 间 的 绿色 
随 接触 角 的 增 大 而 减 小 . 当 接 触角 介 于 0.401 
rad 之 间 时 , 液 滴 厚 度 进一步 增加 , 热 毛 
导 作用 ， 由 于 热 毛细 对 流 的 流动 方向 对答 六 的 BM 
对 流 涡流 的 流动 方向 相反 ， 它 对 妈 膏 满 流 的 产生 具 
有 抑制 作用 ， 因 而 此 时 产生 多 A 和 6 涡流 结构 的 分 
界 Ma 数 随 接触 角 的 增 大 页 增光. 
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Ww 流 结构 
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图 8 三 种 Marangoni 涡流 结构 之 间 的 分 界 Ma 数 
Fig. 8 The Ma number for onset of different kinds of 


Marangoni convection modes versus the contact angle 


4 结 论 


在 考 虚 液 滴 蒸发 冷却 效应 、 气 相 物质 扩散 、 气 液 
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两 相传 热 、 气 液 界面 变形 情况 下 , 数值 模拟 研究 了 固 
定 接触 角 模 型 下 液 滴 莹 发 过 程 中 气 液 两 相 流动 ， 分 
析 了 接触 角 OL Ec BL. Ma 数 对 液 滴 燕 发 过 程 的 影 
响 , 重点 研究 了 液 滴 内 部 Marangoni 对 流失 稳 现 象 。 
得 到 以 下 结论 : 

1) 相同 润 湿 半 径 下 ， 接 触角 0 Be, WOT AR 
寿命 越 短 ， 且 体积 的 2/3 次 方 随时 间 线 性 减少 . 相 
同 接触 角 下 , Bc 数 对 液 滴 蒸 发 寿命 ( 液 滴 体 积 ) 无 明 
显影 响 ; 而 Ma 数 越 大 , 液 滴 无 基 纲 蒸发 寿命 越 长 

2) 当 接 触角 为 0.105~0.175 rad， 随 着 Ma 数 
的 升 高 ， 液 滴 内 部 依次 出 现 热 毛细 涡流 结构 、 单 个 
BM 涡流 结构 及 多 个 BM 涡流 结构 。 当 接触 角 为 
0.175~0.279 rad， 液 滴 内 部 仅 存 在 多 个 BM 涡流 结 
构 。 当 接触 角 介 于 0.279~0.479 rad 之 间 时 , 液 滴 内 
部 存在 藉 涡流 结构 、 多 个 BM 涡流 结构 。 当 接 
触角 0.479 rad 时 , 液 滴 内 部 仅 存 在 热 毛 细 涡 流 


Si 3) 当 接 触角 为 0.105~0.401 rad 时 , 热 毛细 涡流 
结构 与 单个 BM 涡流 结构 、 单 个 BM 涡流 结构 与 多 
个 BM 涡流 结构 、 热 毛细 涡流 结构 与 多 个 BM 涡流 
结构 之 间 的 分 界 Ma 数 随 接触 角 增 大 而 减少 . 当 接 
触角 大 于 0.401 rad 时 , 热 毛 细 涡 流 结构 与 多 个 BM 
涡流 结构 之 间 的 分 界 Ma 数 随 接触 角 增 大 而 增 大 . 
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